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Abstract 

The integrable 3rd-order Korteweg-de Vries (KdV) equation emerges uniquely at linear order in 
the asymptotic expansion for unidirectional shallow water waves. However, at quadratic order, this 
asymptotic expansion produces an entire family of shallow water wave equations that are asymptotically 
equivalent to each other, under a group of nonlinear, nonlocal, normal-form transformations introduced 
by Kodama in combination with the application of the Helmholtz-operator. These Kodama-Helmholtz 
transformations are used to present connections between shallow water waves, the integrable 5th-order 
Korteweg-de Vries equation, and a generalization of the Camassa-Holm (CH) equation that contains 
an additional integrable case. The dispersion relation of the full water wave problem and any equation 
in this family agree to 5th order. The travelling wave solutions of the CH equation are shown to agree 
to 5th order with the exact solution. 

1 Introduction 

We study the irrotational incompressible flow of a shallow layer of inviscid fluid moving under the 
influence of gravity as well as surface tension. Previously Dullin et al. 1111 studied the case without 
surface tension, which in the shallow water approximation leads to the Camassa-Holm equation (CH). 
CH is the following 1+1 quadratically nonlinear equation for unidirectional water waves with fluid 
velocity it (x, t) , 

mt + conix + mrix + 2mux + Tu^xx —O- (1) 

Here m = u — a'^u^x is a momentum variable, partial derivatives are denoted by subscripts, the 
constants and F/cq are squares of length scales and Co = \/gh is the linear wave speed for 
undisturbed water at rest at spatial infinity, where u and m are taken to vanish. The limit 
recovers the KdV equation, |25| . 

Equation ^ was first derived by using asymptotic expansions directly in the Hamiltonian for 
Euler's equations for inviscid incompressible flow in the shallow water regime. It was thereby shown 
to be bi-Hamiltonian and integrable by the inverse scattering transform in the work of Camassa & 
Holm 0. Its periodic solutions were treated by Alber et al (see |2| and references therein). 
Before families of integrable equations similar to were known to be derivable in the general 
context of hereditary symmetries by FOKAS & Fuchsteiner 1161 . However, equation Q was not 
written explicitly, nor was it derived physically as a water wave equation and its solution properties 



were not studied before See |17| for an insightful discussion of how the integrable shallow water 
equation 10 relates to the mathematical theory of hereditary symmetries. 

Equation (0 was recently rederived as a shallow water equation by using asymptotic methods in 
three different approaches by FOKAS & Liu d], by Dullin et al. m] and also by Johnson |^ . 
These three derivations used different variants of the method of asymptotic expansions for shallow 
water waves in the absence of surface tension. We shall derive an entire family of shallow water 
wave equations that are asymptotically equivalent to equation Q at quadratic order in the shallow 
water expansion parameters. This is one order beyond the linear asymptotic expansion for the KdV 
equation. The asymptotically equivalent shallow water wave equations at quadratic order in this family 
are related amongst themselves by a continuous group of nonlocal transformations of variables that 
was first introduced by Kodama |23) . 

Dullin et al. 1111 focused on the integrability of equation Q and its spectral properties. Its 
derivation from Euler's equation in the case without surface tension was briefly described. Here we 
present the necessary details of this calculation. The present derivation also adds surface tension. In 
view of the many papers which have appeared recently on weakly nonlinear shallow water models (i.e. 
I15II22I ') we see our contribution to be the following. By combining a nonlocal Kodama transformation 
with the application of the smoothing Helmholtz operator, we derive several integrable water wave 
equations in the same asymptotically equivalent family. This family includes the CH equation 
the fifth-order Korteweg de Vries equation (KdV5) and the equation of Degasperis & Procesi |S] 
which was recently discovered to be integrable in j^. All these integrable equations are then explicitly 
related to each other, again by means of a Kodama transformation. We clarify the differences among 
the previous derivations of these equations. The equations are discussed with respect to their linear 
dispersion properties. 

In the context of water waves subject to surface tension, interest has recently focused on the KdV5 
equation and its solitary wave solutions. See IIUI for a review. For Bond numbers < a < 1/3 it has 
been shown that these solutions are not true solitary waves which decay to zero at spatial infinity but 
instead they are generalized solitary waves which are characterized by exponentially small ripples on 
their tail. See for example |3] for more explanation. It has been shown by Lombardi I27| that these 
ripples are exponentially small in terms of F — 1 where F — cjco is the Froude number. Numerical 
experiments by Champneys [7| suggest that in the full nonlinear water wave problem there are no real 
solitary waves bifurcating for Bond numbers < a < 1/3. For Bond numbers larger than 1/3, one 
obtains depressions with negative velocity, rather than elevations with positive velocity. 

One may ask whether yet another model equation is needed, if the more rigorous, exact, or numerical 
results described above are already available. Or in more general terms: Which is preferred? An exact 
solution of an approximate model equation, or an approximate solution of an exact equation? Our 
point of view is: If the added cost is small, why not take both? It might be useful to have an equation 
that gives, e.g., more accurate travelling waves than KdV, without the need to go to much more 
elaborate models. Although we shall obtain less information and accuracy than the sophisticated, 
beyond-all-order methods, we shall also pay less. Thus, one can improve the description of the shape 
and speed of the travelling wave without having to resort to these more complicated models. That 
they are still only an approximation to the true solution is, of course, taken for granted. 

Our inclusion of surface tension has a similar motivation. Although the equation 13011 that we shall 
derive has some drawbacks concerning the global properties of its dispersion relation for large fc, it 
still gives improved descriptions for small k and small a . Moreover, the improved solutions are easily 
obtained and analyzed. 

Outline Section |5| recalls the standard dynamics for the shallow water wave elevation following 
Whitham |3U|. We then use an approach based on the Kodama transformation to derive equation 
Q with surface tension in section |3 Section |1] explores the transformations employed to derive wa- 
ter wave equations. We discuss the class of equations which may be related to each other via such 
transformation, and are, hence, asymptotically equivalent. As examples, we discuss the relations of 
equation Q to KdV and other integrable equations. We particularly discuss the relations to KdV5, 
the fifth-order integrable equation in the KdV hierarchy, and to another integrable nonlinear equation 
recently proposed by Degasperis & Procesi [S| and discovered to be integrable in Section |S1 
compares the dispersion relation of Q with that of the full water wave equation. Finally section |S| 
shows that the travelling wave solutions of Q agree to fourth order with the exact travelling waves 
solutions of the full water wave problem. 



2 The 7] equation 



Our derivation of equation proceeds from the physical shallow water system along the lines of 
Whitham I3UI . Consider water of depth h = ho + Tj{x,t) where z = —ho is the flat bottom and ho 
the mean depth, so that 2: = at the free surface in equilibrium. Denote by Uh and the horizontal 
and vertical velocity components, respectively. The z -momentum equation is 

w:th p = a^^^j^, (2) 

where g is the constant of gravity and a is the surface tension. At the free surface, the boundary 
condition is 

Dri 

= Uy where z = rj . (3) 

The irrotational velocity u{x,z,t) — \/(p has horizontal and vertical velocity components Uh = fx 
and Uv — fz ■ The velocity potential ip must satisfy Laplace's equation in the interior. Eq. Q yields 
the kinematic boundary condition at the free surface, 

rjt + fxr)x = fz ■ 

Eq. (|5J can now be integrated to yield the dynamic boundary condition, 

+ ^{<pl + fl) = -gh - -p. 
2 p 

The equations for a fluid are non-dimensionalized by introducing x = Ixx' , z — hoz' , t — {lx/co)t' , 
77 = arj' and ip = {glxCi/ co)f' , where co = y/gho- Being interested in weakly nonlinear, small-amplitude 
waves in a shallow water environment, we introduce the small parameters e = a/ho and 5^ = (ho/lx)^ 
where e > 5^ > e'^ > eS^ > 5^ . When we talk about the linear approximation we lump together the 
terms or order e and , because they are allowed to be of the same order. Similarly, quadratic order 
means quadratic in e and , hence includes terms of order , eS^ , 5* . Contrary to this convention 
the traditional designations 3rd or 5th order KdV refer to the power of 5 , and hence also the maximal 
number of derivatives in the linear terms. Upon omitting the primes and expanding the pressure term 
up to order e^(5^ , the Euler equations and the boundary conditions at the free surface and at the 
bottom are expressed as, 

S'^fxx + fzz = in -1 < z <€T) (4) 

rit + f-fxTix — T^'Pz =0 at z = er] (5) 

V + ft + ^(f^fl + -^fl) - (tS^Vxx = at z = erj (6) 

ipz=0 at z ^ -1 , (7) 

where a = a/{hopCo) is the dimensionless Bond number. The ordering of e and is as specified, 
provided that 0{a) = 1 . 

It is well known that variable transformations of the water wave problem yield the following de- 
coupled equation for the elevation, rj, 

Vt +Vx + 7:'^V'nx + 4<5^(1 — 3a)r]xxx — ^f^V^Vx + <^S^ (777(23 + I5a)r]xrixx 
2d o V 24 

+ -^(5 ~ 3a)r]r]xxx'j + ^^-^i^^ " SOo" - 45a'^)r]xxxxx = 0. (8) 

The derivation up to this order appears, for example, in 1281 . or more recently in |22l without surface 
tension. 



3 Transformation to an integrable equation 

Before embarking on its derivation, we shall discuss the transformation properties of equation Q . First, 
it is reversible, i.e., it is invariant under the discrete transformation u{x, t) —u{x, —t) . Equation (0 
is also Galilean covariant. That is, it keeps its form under transformations to an arbitrarily moving 
reference frame. This includes covariance under transforming to a uniformly moving Galilean frame. 



However, equation Q is not Galilean invariant, even assuming that the momentum, m, Galileo- 
transforms in the same way as velocity, u. In fact, equation Q transforms under 



t ^ t + to , X ^ X + xo + ct , u ^ u + c + uo , m ^ m -\- c + uo , (9) 

to the form 

rut + um^c + 2ux m + (co + uo)mx + 2u^(c + uo) + Tu^xx = . (10) 

Thus, equation (0 is invariant under space and time translations (constants a::o and to), covariant 
under Galilean transforms (constant c), and acquires linear dispersion terms under velocity shifts 
(constant uo ) . The dispersive term uomx introduced by the constant velocity shift uo 7^ breaks the 
reversibility of equation Q. 

Under scaling transformations oi x , t and u , the coefficients of equation Q may be changed. 
However, such scaling leaves the following coefficient ratios invariant, 

CiUxUxx) ■ C{uUxxx) ^ 2 : 1 , (11) 

C{uxxt)C{uux) : C{uUxxx)C{ut) = 3 : 1 , (12) 

where C(/) stands for the coefficient of / in the scaled equation. It is pertinent to mention that the 
above ratios are crucial in the integrability of equation See and 1111 for discussions of this point. 

Equation Q will emerge as being asymptotically equivalent to equation after two steps. First, 
we shall perform a near-identity transformation 

n = n{n) =u + ef[u] + 5^g[u] , (13) 

relating the wave elevation and a "velocity-like" quantity, u. One should not regard u as the original 
fluid velocity, because we will use a non-local term in / which is difficult to interpret in this context. 
Instead, we consider u as an auxiliary quantity in which the equation becomes particularly simple. The 
quantity u agrees at leading order with the fluid velocity and it transforms as a velocity under time 
reversal, spatial reflection and Galilean boosts. To obtain the physically more meaningful quantity r] 
one has to transform back, see below. The functionals / and g in the transformation 11311 are to be 
chosen so that they generate the terms proportional to uux , UxUxx , uuxxx and Uxxx in equation Q. 
Afterwards, we shall apply the Helmholtz operator H = 1 — vS^d^ , which generates the Uxxt term. 
As in 1231 the functional g[u\ is proportional to Uxx and /[it] is a linear combination of and a 
non-local term Uxd~^ , where denotes integration over x. Thus, including the parameter v , there 
are four coefficients in this Kodama- Helmholtz transformation (KH-transformation). These shall be 
chosen so that equation Q emerges, after a rescaling of u , a; and t . 

With these choices, 11311 becomes the Kodama transformation, which depends on three parameters 
ai , 02 and (3, 

rj — rj{u) = u + e(aiu^ + a2Uxd~^u) + S'^ f3uxx ■ (14) 

Terms of degree n in the expansion parameters e and 5^ start contributing at degree n -f 1 in 
the transformed equation. Therefore, no terms of quadratic order in e and 5^ are needed in the 
transformation. Inserting the Kodama transformation 1141 into equation JHJ for the height field rj 
leads to the following terms in asymptotic order. 

0(1) : Ut+Ux 

0{e) : 2aiuut + 2aiuux + a2{uxtd^ u + Uxxd^ u + Uxd^ ut + uux) + -uUx 

0(6^) : fiUxxt + UxxxiP + ^ - (15) 

2\ 9 2,3,2, 0-1 I 2 0-1 \ 3 2 

U(e ) : -aiu + -a2[u + uUxxO u + li^d u) — -u Ux 
z z o 

0<,€5^) ■■ (— + -<T + -(3ai+2a2)(l-3a) + -/3)u,M,, -h 
z4 00 z 



+ 



.1 . „ . , 1 „„ 2^ 



11 3 1 

-a + -(2ai + 3a2)(l — 3o-) + -f3)uUxxx + ^^2(1 — 'i(T)uxxxxd~ i 



0{5*) : (-/32(l-3a)-f — (19-30cr-45cr''))u,,,,, . 
o 360 



As before, we expand the time derivatives to linear order as 



3 1 

Ut = —Ux — -eUUx — (1 " 3cr)Ua;a:x , (16) 

3 2 3 1 2/, o \ 

Uxt = —Uxx — ~ 2*"""^^ ~ 6 ■io-jUxxxx , 

9 _ 3 , . 

xxt — Uxxx ^^UxUxx ^^UUxxx ■ \ '/ 



This expansion generates higher order terms, leading to 

Oil) : ut+Ux 
3 

0{e) ■ -uux 
^2^ 1 



0(5") : -{l-io)uxxx (18) 



0(6^) : (la, + \a2-l)u^ux (19) 

0(65^) : AUxUxx + BuUxxx (20) 

0(<5'') : ^(19-30a-45o-')?i,,,., . (21) 



where in 1201 we defined 



- 23 5 1 - '5 1 1 

^ = — + -(7+ -(2ai +a2)(l - 3cr) - 3/3 and B = cr + -02(1 - 3cr) . 

24 8 2^ ' ' 12 4 2 ^ ' 

The first step of the derivation is now complete. In the second step, applying the Helmholtz operator 
H = 1 — tyS^dx introduces the coefficient ly and creates terms with two more x derivatives. However, 
the terms of order 0{e^) are left unchanged. These terms are proportional to u^Ux and they must 
vanish for equation Q to emerge. Thus, the application of the Helmholtz operator restores the needed 
Uxxt term that had previously been eliminated. Alternatively, the same equation (0 can be obtained 
by splitting the time derivative, that is, by partially substituting the time derivative Uxxt in 1151 by 
its asymptotic approximations 117L along the lines of Benjamin et al. |2] in their study of the BBM 
equation. 

The order O(e^) coefficient in expression 119t will vanish, provided the parameters ai and 02 are 
chosen to satisfy 

4ai + 2a2 = 1 . (22) 

The order 0{5^) terms receive an additional contribution that arises from applying the Helmholtz 
operator H = 1 — uS^dx to the terms of order 0{S^) . The entire combination at order 0(5'') must 
vanish, for the final equation to possess no fifth-order derivative term, Uxxxxx ■ This requirement 
determines v as 

1 19 - 30cr - 45(j2 

v = . (23) 

60 1-3(7 ^ ' 

In what follows, we shall consider the coefficient v to be given by this function of surface tension, 

a. Note: this removal of the highest order term was made possible by introducing the additional 

parameter u via the Helmholtz operator. The remaining terms containing free parameters 02 and /3 

are of order eS^ and they combine additively as 

- 9 - 3 

(A — -u)UxUxx + {B — -v)uUxxx ■ 

To ensure equivalence to Q except for scaling we need the relative coefficients to appear in the ratio 
ljll|l . so that 

(i„ 9^) 1^)^2:1. (24) 

In addition we also need to satisfy 11211 . so that 

2 ^ 2 ' 

These two conditions imply B = 2u and A — llv/2. As a result we finally obtain the equation 

Ut — fS^Uxxt + Ux + ^euux — ^eS'^u{uUxxx + 2UxUxx) + ^^(^ ^ ^ ^ ^cr)uxxx = , (25) 



which can be rewritten in terms of m = u — v5 Uxx as 

mt + mix + ^{urrix + 2mUx) + <5^(- — -(j)uxxx = . (26) 

Thus, the coefflcients in the Kodama transformation 114^ that yield this equation are, 

7 1 2-3cr _ 

^ 20^^5(1^3^' (2^) 

1 2 2-3a 

= -5+^5(13^' (2^) 

^ 30 30 1 - 3cr ^ ' 

Returning to physical variables, where u and (fix have units of ga/co — coa/ho, followed by an 
additional scaling of m — > 2ii, gives equation 1261 the canonical CH form, 



mt + corux + umx + 2mux + Fuxxx = , (30) 

in which m = u — vh^Uxx and F = co/i^(l — 3ct)/6. The parameters o? and F in CH are now 
understood in terms of physical variables as 

2 ,2 ^2 1 19 - 30a - 45a' ^ coh? 
"="'^='^60 1^3^ ' (l-3a). (31) 

The parameter F changes sign when the Bond number a crosses the critical value, 1/3. For later 
reference, we record the value of a > for which o? vanishes, as 

ac = -i + :j|\/30 ~ 0.39696 > \ . (32) 

In the special case cq = F = equation 13011 is called "the peakon equation." This equation 
supports peakons as solitary wave solutions whose derivative is discontinuous at the extremum. These 
solutions were introduced and discussed by Camassa & Holm [5]. The peakon equation has many 
exceptional mathematical properties that arise from its interpretation as geodesic motion in the Euler- 
Poincare variational framework, as explained in I20| . However, the peakon equation cannot be derived 
as a water wave equation in a weakly nonlinear shallow approximation from the Euler equation by the 
present technique. This is because neither a Galilean transformation, nor an appropriate splitting can 
eliminate both of the linear dispersive terms in equation 13011 simultaneously. One is always left with 
a residual linear dispersion. 

Johnson 1221 has recently derived the dispersive CH equation 13011 as a shallow water wave equation 
by using the same asymptotic expansion. However, two key steps in the asymptotic derivation are 
addressed in quite different ways in 1221 and here. Firstly, the ratios of the coefficients of the equation 
must be adjusted to ensure integrability. In 1221 . this is achieved by using the height at which the 
velocity potential is evaluated as a free parameter. Instead of the height, we use the free parameters 
in the Kodama transformation I14II to obtain the desired ratios. 

Secondly, the fifth-derivative term Uxxxxx of order 0(S'^) must be removed. The present approach 
allows the value of the free parameter f in the Helmholtz operator to be chosen to cancel out this 
term. In 1221 . this term is simply omitted. In Section 4 we shall use our approach to show that the CH 
equation is asymptotically equivalent to the KdV5 equation, which involves this fifth-order derivative. 

In order to compare predictions in terms of physically measurable quantities, the solutions for the 
velocity-like variable u must be transformed back to the elevation field rj by using 1141 . However, the 
derivation not only used the transformation 1141 . it also involved applying the Helmholtz operator. 
Therefore, one should verify that it is sufficient to simply invert 1141 . Fortunately, when the inverse 
transformation u — uijf) of the same form as 1141 with u and rj interchanged is substituted into 
the final equation 1261 . we find that the coefficients just reverse their signs. We conclude that (0 is 
equivalent to the shallow water wave equation Q up to and including terms of order 0{5*) . 



4 The Kodama-Helmholtz transformation 

In the previous section, we transformed the shallow water wave equation JSJ into the CH-equation Q 
by means of a Kodama transformation 1141 and an application of the Helmholtz-operator. Now we 
describe the class of equations that can be derived from the shallow water wave equation JHJl by any 



such sequence of transformations, which we will refer to as KH-transformation. This class contains 
(at least) four integrable equations, one of which is the CH-equation. Moreover, all equations in this 
class are not only related to l|SJ; they are also related to each other by such transformations and are 
therefore asymptotically equivalent. 

Similar transformations have been used by FOKAS & Liu |14) that included an additional term of 
the form xut , but no Helmholtz-operator. In this case, the class of equations even includes the (3rd 
order) KdV equation. Unfortunately, though, this term is not uniformly bounded, so we shall decline 
to use it. Its unboundedness is a problem when transforming travelling wave solutions which move 
towards x = ±00 . Moreover, use of the term xut changes the dispersion relation. 

4.1 Invariance of the dispersion relation 

The Kodama transformation l|14|l does not change the dispersion relation. To see this, one may observe 
that only terms linear in rj or its derivatives produce linear terms in the transformed equation. Simi- 
larly, applying 1141 to nonlinear terms in an equation produces only nonlinear terms in the transformed 
equation. Therefore, in proving invariance of the linear dispersion relation under H14|l . we may restrict 
to a transformation 77 — u + eL{u) in which L is a linear differential operator with constant coefHcients 
and the linear equation to be transformed is rjt — M{r^) . To first order, we then have Ut = M{u) and 
the full transformation gives 

Ut + eL{ut) = M{u + eL{u)) . 

Now the first order equation may be used to eliminate the time derivatives which are not of order zero, 
thereby yielding 

Ut + eL{M{u)) = M{u) + eM{L{u)) . 

If M and L commute, as they always do when they have constant coefficients, the final answer is 
ut = M{u) . Consequently, the Kodama transformation 11411 leaves a linear equation unchanged. Note 
that including the xut term in the transformation would in general cause the operators to no longer 
commute, so that the dispersion relation would be changed, as previously claimed. 

The second step of the transformation is the application of the Helmholtz operator H — 1 — i^S^d^ . 
As we have just seen, the Kodama transformation leaves the linear part of the equation unchanged. 
Applying the Helmholtz operator to an equation does change the linear part, but it still leaves the 
dispersion relation unchanged. To see this, let the linear part of the equation be given by Ut — M{u) . 
The new equation is H{ut) — H[M{u)). If H and M are linear with constant coefficients this gives 
H{u)t — M{H{u)) so that with the definition m — H{u) we obtain mt — M{m) , which has the same 
dispersion relation. This does not hold, however, if we truncate higher order terms in H{M{u)). If 
we truncate, then the dispersion relation will agree up to the order of truncation. For example, the 
dispersion relation for (0 is a rational function, which differs from the polynomial dispersion relation 
obtained from JHJ. However, by the above argument the two dispersion relations agree up to the desired 
order. 

4.2 Range of the KH-transformation 

We now investigate which equations may be transformed into each other by a KH-transformation. The 
class of equations that (under some additional conditions to be derived) may be transformed into each 
other by KH-transformations is given by 

F{ut,Ux,Uxxx,Uxxt,Uxxxxx,UUx) + G{u^Uj:,UUxxx,UxUxx) =0, (33) 

where F and G are linear and the coefficients of Ut , Ux , uux , Uxxx are nonzero. This means that 
each equation has a KdV-kernel. In addition, we assume that the terms are ordered as in the previous 
section: Every u has weight e and every x (or t) derivative has weight 5. However, for simplicity, 
we do not explicitly display the weights in the following, even though they are used for truncation at 
the usual order. The general procedure consists of two parts. In the first part the coefficients in F are 
normalised in three steps. In the second part the Kodama transformation is used to adjust the terms 
in G . The details of the calculation are similar to those of the previous section and are therefore not 
given. 

Equations are considered to be equivalent when they differ only by a transformation of the form 



{t, X, u) — > {rt, ^a; — ftt, vu) . 



(34) 



Scaling of u, t, and the equation allow one to set the coefficients of Ux , ut, and uux to arbitrary 
non-zero values. After such a scaling, the equation takes the form, 

Ut + Ux + UUx + fsUxxx + fAUxxt + fsUxxxxx + QlU^Ux + g2UUxxx + Q^UxUxx = . (35) 

The next step is to remove the Uxxt term by eliminating the t derivatives using the equation itself. 
This operation produces an equation of the form, 

Ut+Ux-^ UUx + (/3 — f4,)Uxxx + (/s — fi{h ~ fi))Uxxxxx + 

giU^Ux + (32 — fA)uUxxx + {g-J, — 3f4)UxUxx = . (36) 

In the third step, the coefficients of Uxxx and Uxxxxx are normalised using 13411 again, while keeping 
the coefficients of itt , Ux , uUx all fixed at unity. First ^ is used to normalize the ratio between the 
coefficients of Uxxx and Uxxxxx to ±1, assuming that both are non-zero. The relative sign of these 
coefficients cannot be changed. Dividing the equation by r allows both coefficients to be set to unit 
magnitude. In order to keep the other coefficients in the linear function F equal to unity, the special 
transformation of the form H34|l reads {t,x,u) — > {Tt,£^x -\- [t — ^)t, £,u/t) . The above three steps 
(transformation I34I I to find 13511 . elimination of Uxxt to find 13611 and again transformation 13411 to 
find 1371 ') can be performed for any equation in the class 13311 . Therefore, we shall only be concerned 
with the equivalence under Kodama transformation of equations of the form 

Ut +Ux + feUUx + fsUxxx + fsUxxxxx + giU^ Ux + g2UUxxx + (jsUxUxx — , (37) 

where — 1, /s = 1, or 0, and /s = ±1, or 0. fi~0 appears when /a — fi in 13511 . The parameter 
/e is kept in the notation to trace the influence of the uUx term. The terms in the linear function 
F are unchanged by the Kodama transformation because such linear terms are unchanged by a near 
identity transformation and uux is the lowest order nonlinearity. Therefore two equations can only 
be equivalent if their normalised coefficient /s agrees. Recall that Uxxt has already been eliminated, 
so that if, e.g., the above steps are applied to the CH-equation, the resulting equation does possess a 
non-zero fifth derivative term. 

The terms in the linear function G in 1331 can be adjusted by the Kodama transformation. This 
generates coefficients, 

gi = 51 + /6(2ai +«2)/2, 
32 = 52 + 3a2/3 , 
53 = 53 + 3/3(201-^02) -2/3/6 

The cji are linear in the coefficients ai , 02 , and jS of the Kodama transformation. A sufficient 
condition for the solvability of this linear system is /3/g 7^ . Hence, if the coefficients of Uxxx and uUx 
in 1371 are non-zero, then any value of gt can be achieved by some Kodama transformation. Since the 
coefficient of uux in 1331 is assumed to be non-zero the only additional condition is /3 = /3 — /4 7^ 
in 1351 . Note that this condition was already necessary in order to obtain 1371 . For example, any 
equation with positive relative sign of /3 and /s in 1371 is equivalent to the water wave equation © 
with (T = 0. In particular, this shows that the integrable KdV5 can be obtained, see 1261 . Similarly, 
different higher-order, but non-integrable, extensions of the KdV equation can also be obtained, e.g., 
those introduced by Champneys & Groves Notice that the relative sign of /s and /s in (l?7l is 
negative for 1/3 < ct < (Tc, see 1321 . The condition /a = appears for a = 1/3, while /s = for 
(J = cTa ~ 0.39696. 

To arrive at equations with a Uxxt term, the Helmholtz operator is applied. The Helmholtz operator 
fits nicely into the above procedure, because it is the inverse of the elimination of the Uxxt term, at the 
order considered. More precisely, if the Helmholtz operator 1 — fid^ is applied to 1361 then equation 
1351 is recovered. Therefore, each of the three steps leading to the normalised form 1371 can also be 
inverted, and any two equations with equal /s are asymptotically equivalent. The CH-equation has 
/s — sign(a^r) , even though /s = originally. Moreover, /a = sign(r) , so that the relative sign is 
that of o? . From these equations, it follows that the vanishing of /a implies the vanishing of /s . This 
is not true in the water wave equation. That is, o — 1/3 does not make the coefficient of Uxxxxx 
vanish. This explains why the CH-equation is not a good model for ct = 1/3. 

In general, all the terms in the linear function G in equation 1331 may be removed by a Kodama 
transformation. However, it is not possible in general at the same time to remove the fifth order 
derivative. One may, however, trade the fifth order derivative for the Uxxt term. Therefore, possibly 
the simplest representative equation in the class 1331 has the form 

Ut +Ux + UUx + Uxxx + K.Uxxt — , 



(38) 



where small positive/negative n gives the different signs for /s . The case with positive k, results from 
the BBM equation of Benjamin et al. by a Galilean transformation and velocity shift (which is 
not included in the transformation group 1341 '). It is known that the BBM equation is not integrable. 
Thus, the integrable CH equation arises as a Kodama-Helmholtz near-identity transformation of the 
non-integrable BBM equation, after a Galilean transformation and velocity shift. 

4.3 Examples 

4.3.1 Asymptotic equivalence of CH and KdV5 equation 

The previous section constructs a transformation from the water wave equation to the CH equation. 
The general argument of this section shows that it is also possible to transform the CH equation, for 
example, into the KdV5 equation. To this end, we first expand the time derivative in the u^xt -term 
using the equation itself and then apply a transformation of the form 

u = V + e{aiv^ + a2Vxd^^v) + 5^ (3vxx 

Choosing the values in this Kodama transformation as 

ai = ^, a2 = 2^, /3 = 2a^ (40) 
transforms the CH equation Q into the integrable KdV5 equation 

Vt + CoVx + ivVx + 5 {vVxxx + 2VxVxx) C? + " ^ ""^ + V [aVxxxxx + Vxxx) = . (41) 

We conclude that Q is asymptotically equivalent to the integrable KdV5 equation, and both of 
them are equivalent to at order 0{S*) . However, the equivalence of Q to the KdV5 equation 
breaks down in the limit F — > 0, because both the transformation and the resulting equation are 
singular in the limit F ^ 0. Therefore, the peakon equation cannot be transformed into KdV5. 

Using the additional parameter v supplied by the Helmholtz operator allows for the removal of the 
highest order term while preserving the dispersion relation, which is unchanged by applying a linear 
operator to the equation. One advantage of the CH equation over the asymptotically equivalent KdV5 
equation is that it is numerically easier to integrate because it does not contain the fifth derivative. 
This is in accordance with the general smoothing effect of the Helmholtz operator. 

Different variants of higher order KdV equations have been derived in the literature. In this context, 
we mention who obtained a 5th order KdV equation by expanding the variational formulation of the 
water wave problem. Again, their equation is asymptotically equivalent to CH by a certain Kodama 
transformation. From our point of view, however, the virtue of the power of the Kodama-transformation 
is that it allows one to select integrable equations as the model equation. The equation derived in |B] 
is not known to be integrable. 



(39) 



4.3.2 Asymptotic equivalence to the 6-equation 

Recently a new variant of Q has been introduced by Degasperis & Procesi [H| as 

mt + urux + bux-m = cqUx — Fuxxx , (42) 

where b is an arbitrary parameter. The cases b — 2 and 6 = 3 are special values for this equation. 
The case 6 = 2 restricts it to the integrable CH equation. The case 6 = 3 is the DP equation of 
Degasperis & Procesi [H], which was shown to be integrable in The two cases CH and DP 
exhaust the integrable candidates for (1421 . as may be shown using either Painleve analysis, as in 
or the asymptotic integrability test, as in The b-family of equations (I42II was also shown in 1291 to 
admit the symmetry conditions necessary for integrability only in the cases 6 = 2 for CH and 6 = 3 
for DP. 

We shall show here that the new integrable DP equation can also be obtained from the shallow 
water elevation equation (|HJ by an appropriate Kodama transformation. The derivation in the previous 
section is essentially unchanged up to equation 1241 . The two scaling relations 1111121 now read 

(i_9^):(B_3^) = 6:1, 

^u:{B-^u) = 6 + 1:1. 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 hk 0.2 0.4 0.6 0.6 1 1.2 1.4 1.6 hk 




0-9 0.2 0.4 0.6 0.6 1 1.2 1.4 1.6 hk o 0.2 0.4 0.6 0.6 1 1.2 1.4 1.6 hk 



Figure 1: Normalised dispersion relations for a = 0.0,0.25,0.35,0.8 as a function of hk. The full line is 
the exact phase speed the dashed line is the approximation (|^ and the short dashed line is that of 
KdV5. 



These two conditions imply 



5 36 + 2 ^ . 3 46 + 3 
B = V — - and A = v — - 



26+1 2 6 + 1 

The resulting Kodama transformation of the form 1141 with coefficients a'l , ct^ , and (3' are 



Qi = Ql + 3A 
a'2 = 02 — 6A 
13' = /9 - (1 - 3o-)A where 
_ 6- 2 45cr^ + 30cr - 19 
6 + 1 360 ' 

Therefore any 6 7^ —1 can be achieved by an appropriate Kodama transformation. Note that when 
cr = (Tq , see 1321 , then = , hence A = is independent of 6 . After tliis transformation 1421 is 
obtained by further scaling the new dependent variable u by the factor 6 + 1 . See 1211 for discussions 
of equation 14211 in which 6 is treated as a bifurcation parameter when co = and F = . 

The value 6 = —1 is excluded, not from a deficiency of the KH-transformation, but because the 
term uu^ is not removable in the linear function F in equation 13311 . As we have shown above, the 
KH-transformation does not affect this term, and so it cannot be removed from the 77 equation 
That is, the case 6 = — 1 is not within the range of the KH-transformation. 

We conclude that the detailed values of the coefficients of the asymptotic analysis at quadratic order 
hold only modulo the KH-transformations, and these transformations may be used to advance the anal- 
ysis and thereby gain insight. Thus, the KH-transformations may provide an answer to the perennial 
question "Why are integrable equations found so often, when one uses asymptotics in modelling?" 



5 Dispersion relation 



The interplay between the local and nonlocal linear dispersion in the CH equation 130|l is evident in 
the relation for its phase velocity, 

for waves with frequency u and wave number k linearized around u = Q. For F < , short waves and 
long waves both travel in the same direction. Long waves travel faster than short ones (as required 
in shallow water), provided F < 0. Then the phase velocity lies in oj/fc G (co — F/a^,co]. At low 
wave numbers, the constant dispersion parameters c? and F both perform rather similar functions. 
At high wave numbers, however, the parameter properly keeps the phase velocity of the wave from 
becoming unbounded, and the dispersion relation is similar to the original dispersion relation for water 
waves, provided that the surface tension vanishes o" = 0. The remarkably accurate linear dispersion 
properties close to « give the CH equation a clear advantage over the KdV equation (provided a 
is not close to 1/3, see Fig. 0. 

For the peakon equation, cq = F = and linear dispersion is absent. For nonvanishing surface 
tension, the true dispersion relation for shallow water waves is unbounded for large wave numbers, 
whereas the dispersion relation of equation saturates to the asymptotic value co — F/q . 

Consequently, (0 is inferior to KdV5 for non-zero surface tension, with respect to traveling wave 
solutions exhibiting small tail waves. Unboundedness of the linear dispersion relation of KdV5 for high 
wave numbers allows resonances to occur between a supercritical solitary wave and high-wavenumber 
linear waves. These resonances give rise to exponentially small ripples at the tails of the solitary wave, 
in accord with the water wave solution of the full Euler equation. See, for example, 111 1191 ITHl and I27| 
for more discussion and analysis of these fascinating resonances. Nonetheless, traveling wave solutions 
of CH and KdV5 are asymptotically equivalent, and they both agree asymptotically with traveling 
wave solutions of the full water wave problem, see sec. |S] 

The connections to the physical parameters a = oi{a) and F = F((j) are defined in equations 
1311 1. The asymptotic equation (I30II was derived from the water wave equation by means of the KH- 
transformation. As explained in sec. 14.11 the dispersion relation (I43II matches the dispersion relation 
for water waves up to quintic order. In comparison, the dispersion relation for water waves, when 
developed for small wave number k yields 



w /l + ah^k^ , , , ,,,, 

— = \ , , tanh/ifc (44) 

cok y hk 

^ 1 - i(l - 3cr)/i^fc2 + ^(19- 30a -45cr^)/i''fc* (45) 

^ 1 - i(l - 3cr)/i^fc^(l - i/ft^/fc^) . (46) 
6 

Therefore, the dispersion relations are in agreement up to 5th order in dimensionless wavenumber hk . 
For the particular value a ~ 0.069, the dispersion relations agree up to 7th order. 



6 Travelling wave solutions 

One may also ask how the solutions of equation (0 compare to the usual KdV solitons and to real 
water solitary waves. An expansion of the form of the traveling wave in the full Euler equations is 
given by Grimshaw ^S], and to higher order by Fenton ^ The result given in [TS] is normalised 
so that the highest point of the wave at s = is unity. 

By direct substitution of a series in sech^ into 1251 1. one finds 

it(s) = asech^bs + ^ea'^sech'^fos , 

where = Sae/{45^) and c = 1 + ea/2 + i|e^a^ in s = x — ct. Applying the (inverse) Kodama 
transformation gives 

rj{s) — a{l + — ea)sech^6s + -ea^sech*6s . 

For simplicity, we restrict the calculation to the case without surface tension, i.e., a = Q. After 
normalising the height at the crest to unity, we find perfect agreement up to order eo? with the exact 



^For only moderately steep waves, and for the corresponding cnoidal waves, Fenton has shown that the 1972 approach 
gives a very poor approximation to the velocity field, see Fenton (1990) [H. Peregrine (private communication)]. 



result by Grimshaw |18| . This shows that the traveUing wave solutions of Q agree with those of the 
exact Euler equations up to the order we are considering, which is the optimal result. In particular, 
the traveling waves are narrower and slower (in this normalisation) than the KdV soliton, in agreement 
with experimental findings of Weidman & Maxworthy The solution also agrees with the one 
found in 1261 by solving the 5th order KdV equation. 
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